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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 



1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 



(fill out in typewriter or pencil, do not use ink) 

Program No. Date_ 

Program Name: 



1. Does the abstract adequately describe what the program is and what Yes No. 

it does? 

Comment 

2. Does the program do what the abstract says ? Yes No_ 

C ommen t 

3. Is the Description clear, understandable, and adequate? Yes No 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No_ 

Comment 

Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable? Yes No_ 

Comment 

5. Does the source program compile satisfactorily (if applicable)? Yes No_ 

Comment 

6. Does the object program run satisfactorily? Yes No_ 

C omment 

7. Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

C omment 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No_ 

Comment 

9. Were all necessary parts of the program received? Yes No_ 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt Company 
Data Corporation 

7500 Old Xenia Pike Address 
Dayton, Ohio 45432 



User Group Code_ 



THIS REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVITED 
TO PARTICIPATE IN THIS EVALUATION. , 

' 11/09/64 
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Subroutine Crout 
for the solution of simultaneous linear equations 




Joyce Fodor 

Engineering Computing Laboratory- 
University of Wisconsin 
Madison, Wisconsin 53706 

User: 31 55 
Dates 12-22-6J4 
Program : 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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Program Abstract 

Title ; Subroutine Crout for solution of simultaneous 

linear equations 

Author ; Mrs. Joyce Fodor 

Engineering Computing Laboratory 
University of Wisconsin 
Madison, Wisconsin 53706 

Date : 12-22-6k 

Users Group Code ; 31 55 

Direct Inquiries to : Prof. C. H. Davidson 
Director, 

Engineering Computing Laboratory 
University of Wisconsin 
Madison, Wisconsin 53706 
Phone 608-262-3892 

Descr ip tion/ Purpose 

This subroutine will solve a system of simultaneous linear equations 
using the Crout method. The number of equations is limited only by the 
size of the main program and the core available. 

Specifications 

1 . 194 1 620 able to use FORTRAN II 

Programming type FORTRAN II subroutine 

Language used in writeup English 
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Deck Labelling Sheet 

Deck No. Name. Card numbers 

1 Subroutine Crout 1 — 60 
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Subroutine Croat for solution of simultaneous linear equations 

Mrs. Joyce Fodor 

Engineering Computing Laboratory 

University of Wisconsin 

Madison, Wisconsin 53706 

User Code 31 55 

Program Developed 0ns 

1 . n 1 620 MODEL II - 60 K 

2. Indirect addressing 

3. 1 disk drive 

it. TNS, TNF, and MF special instructions 

5. Floating point hardware ,^ - ro» :•.•'* \ 

Machine Configuration Required; 

1 . im 1 620 capable of using FORTRAN II 

System used: 

FORTRAN II 



Program Description 

Given the set of equations [A] [X] = [B] the subroutine will solve 
them using a modified Crout method. . In reducing the coefficient deter- 
minate the subroutine searches out the largest coefficient on or below 
the major diagonal in the column being reduced. It then interchanges 
the rows, placing this coefficient on the diagonal before continuing. 
This reduces the round-off error in the calculation. The main program 
must read in or calculate the coefficients of the equations. 

When writing the main program the size of the dimensioned arrays 
of the coefficient matrix (A in the subroutine), the constant terms (B), 
and the solution (X) must agree between the main program and the sub- 
routine. The three items above must also be put in common in that order. 

To call the subroutine a statement 

CALL CROUT (N1, V1 ) 

should be used where N1 is a fixed point variable or constant equal to 
the number of simultaneous linear equations and V1 is a floating point 
variable name. 



* F. 5. HiWebrand, Introduction to Numerical Analysis, McGra»/-Hill, 
New York, 1?56, pp. U29-WJ,'. 
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When returning from the subroutine V1 should be checked. If 
VI = 0. a solution has been reached. If V1 » 1 . , the subroutine 
has found that the set of equations is probably dependent and a 
solution has not been reached. This indicates that the largest term 
in any column on or below the major diagonal is less than 1.X10-30. 
Unless the coefficients of the matrix are very small this condition 
would indicate that the equations are dependent. Therefore, if the 
coefficients; of the matrix are very small to begin with, the constant 
in statement 0021+01 lines should be made smaller. 

The operating instructions are the same as for any FORTRAN II 
subroutine. There are no sense switches used. 
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SUBROUTINE CROUT(N.V) 


1 










DIMENSION A(9,9) »B(9).X(9) 


2 




DIMENSION A(9»9)»B(9)»X(9). 2(9.9) 


1 




COMMON A.B.X 


3 




COMMON A.B.X 


2 




NM1=N-1 


4 


18 


READ100.N 


3 




DO 2 J = l »N 


. 5 




IF ( N ) 106 » 106 .107 


4 




JP1= J+l 


6 


106 


PUNCH200 


5 




JM1=J-1 


7 




STOP 


6 




DO 6 I=J»N 


8 


200 


FORMAT ( / 2X27HTH I S IS THE END OF THE DATA) 


7 




ASUM=0 . 


9 


100 


FORMAT (12) 


8 




IF (JMU6.6.7 


10 


107 


CONTINUE 


9 


1 


DO 9 K=1»JM1 


11 




READ 101. ( (A(I.J).J=1»N).I=1.N) 


10 


9 


ASUM= ASUM+A ( I »K ) *A < K » J ) 


12 


101 


FORMAT (5E15.0) 


11 


6 


A t I » J ) =A( I » J)-ASUM 


13 




READ 102. <B( I ) »I=1»N) 


12 




AMAX=A(J»J) 


14 


102 


FORMAT (5E15.0) 


13 




IMAX=J 


15 




DO 108 I=1»N 


14 




IF ( JP1-NJ20. 20.21 


16 




DO 108 J=1»N 


15 


20 


CONTINUE 


17 


108 


Z ( I . J ) =A ( I .J) 


16 




DO 1 I=JP1»N 


18 




CALL CROUT (N.V) 


17 




I F ( ABSF ( AMAX ) -ABSF ( A< I.J) )) 3.1.1 


19 




IF (V) 120.120.10 


18 


3 


AMAX=A (I.J) 


20 


120 


CONT INUE 


19 




IMAX=I 


21 




DO 103 1=1 »N 


20 


1 


CONTINUE 


22 




PUNCH 104. 1 .X( I ) 


21 


21 


CONTINUE 


23 


104 


FORMAT (12H X(I2.5H)= E15.5) 


22 




IF ( ABSF ( AMAX ) -1 .E-30 ) 10. 10 .4 


24 


103 


CONTINUE 


23 


10 


V = l» 


25 




DO 109 1 = 1. N 


24 




RETURN 


26 




SUM=0. 


25 


4 


DO 5 K = 1»N 


27 




DO 110 J=1»N 


26 




ASAVE=A( IMAX.K) 


28 


110 


SUM=SUM+2 ( I . J ) *X ( J ) 


27 




A( IMAX.O = A( J.K) 


29 


109 


PUNCH lll.I.SUM 


28 


5 


A( J.K) =ASAVE 


30 


111 


FORMAT ( 3X9HCONSTANT ( I 2 . 5H ) = E15.5) 


29 




ASAVE=B ( IMAX ) 


31 




GO TO 18 


30 




B( IMAX ) =B ( J ) 


32 


10 


PUNCH 105 


31 




B( J)=ASAVE 


33 


105 


FORMAT (/2X22HTHE MATRIX IS SINGULAR/) 


32 




1 1 = J 


34 




GO TO 18 


33 




J1=JP1 


35 




END 


34 




IF (JP1-NJ22. 22.23 


36 








22 


CONTINUE 


37 










DO 8 J2=JP1.N 


38 










ASUM=0. 


39 










IF (JM 1)8. 8. 11 


40 








11 


DO 12 K=1.JM1 


41 








12 


ASUM=ASUM+A< 1 1 . K ) *A ( K > J2 ) 


42 








8 


A(I1»J2) = (A(I1»J2) -ASUM ) / A ( 1 1 . 1 1 ) 


43 








23 


CONTINUE 


44 










ASUM=0 . 


45 










IF (JMD2.2.13 


46 








13 


DO 14 K=1.JM1 


47 








14 


ASUM= ASUM+A ( 1 1 . K )*B( K ) 


48 








2 


B ( 1 1 ) = (Bill) -ASUM ) / A ( 1 1 . 1 1 ) 


49 










DO 15 J=1.N 


50 










I 1=N-J+1 


51 










1=11+1 


52 










ASUM=0. 


53 










IF ( Il-N117.15.15 


54 








17 


DO 16 K=I»N 


55 








16 


ASUM = ASUM+A( I1.K)*X(K) 


56 








15 


X( I1)=B( 11 J-ASUM 


57 










V = 0. 


58 










RETURN 


59 










END 


60 
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ft****************** EQUATIONS ****•> 



.4QX1 
, 64X 1 
.6 XI 



i 67X2 
, 67X2 
, 33X2 



.6X3 
.8X3 
0.X3 



= 0. 

=4000. 
= 0. 



##*###*#•»•#***■*•*#*■***** j NPUT ****** ****************** 



3 

-.48 



0. 



-.67 
.6 

4000. 



. 6 

-.33 



******************** KESULT5 **********************■! 



X ( 
X( 
X( 

CONSTANT ( 
CONSTANT ( 
CONSTANT ( 



1 ) = 

2 ) = 
3) = 
1 ) = 

2) = 

3 ) = 



970 .30300E+00 
176.41870E+01 
274.62511E+01 
000.0OOO0E-99 
399.99999E+01 
900.00000E-07 



a 



